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Abstract 

Recent results show that the peak transient response of a linear system to 
bounded energy inputs can be computed using the energy-to-peak gain of the 
system. However, analytically computed peak response bound can be 
conservative for a class of bounded energy signals, specifically pulse trains 
generated from jet firings encountered in space vehicles. In this paper, shaping 
filters are proposed as a methodology to reduce the conservatism of peak 
response analytic bounds. This methodology was applied to a realistic Space 
Station assembly operation subject to jet firings. The results indicate that 
shaping filters indeed reduce the predicted peak response bounds. 


1. Introduction 

Assembly of the International Space Station requires robotic manipulation of various 
payloads. The payloads will be manipulated by the Shuttle Remote Manipulator System 
(SRMS), the Space Station Remote Manipulator System (SSRMS). During assembly operations, 
it may be necessary to execute corrective attitude control thruster firings. The operational plan 
requires stopping the robotic payload operation and activating the SRMS/SSRMS joint brakes 
before firing attitude control thrusters. Hence, to evaluate the feasibility of planned robotic 
operations, it is necessary to show that the joint brakes do not slip during thruster firings. 
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To evaluate the feasibility of thruster firings when the robotic systems are in the brakes- 
on mode requires evaluating the maximum joint torque and comparing that with a brake slip 
limit torque. Typical attitude control firings include bipolar pulse trains representing attitude 
maneuvers and unipolar pulse trains for attitude hold or momentum desaturation. The standard 
approach to evaluate the peak brakes-on joint torque is via time-domain simulations. This 
approach can be time consuming and inefficient since time-domain simulations can only 
calculate the peak response to one pulse train at a time. Thus, it is difficult to map the time- 
domain simulation results to a family or class of pulse trains. 

Recent results [1,2] in linear operator norm theory show that the peak transient response 
of a linear system to bounded energy inputs can be computed using the energy-to-peak gain of 
the system. This computation requires the solution of a Lyapunov equation of order equal to the 
order of the system. The energy-to-peak gain theory also provides the worst-case input that will 
result in the transient response reaching its peak. It is noted that the bounds provided by this 
theory are tight, that is, the bounds are achievable by a certain input excitation. 

The class of all energy-bounded input signals considered in the energy-to-peak gain 
theory is a very broad one for practical applications. Hence, the analytically computed peak 
response bounds can be conservative for practical classes of bounded energy signals, such as, 
pulse trains generated from jet firings encountered in space vehicles. In this paper, shaping 
filters are proposed as a methodology to reduce the conservatism of peak response analytic 
bounds computed by the energy-to-peak gain theory. A filter design procedure corresponding to 
a given class of input signals is presented for conservatism reduction, and the methodology is 
illustrated using two types of input signals for a realistic example involving Space Station robotic 
assembly. 

The notation to be used in this paper is as follows: ||u||^ denotes the energy (or 2-norm) of 
a signal u, that is. 


IML = (J «(O r u(t)dt) 112 

o 

and |w|| denotes the peak (or infinity) norm of a signal u, that is 

||u||, = max |u,(r)| • 

II H L„ j 

For a matrix A, IIAII denotes the maximum singular value (or induced 2-norm) of the matrix. 

2. Worst-Case Energy-to-Peak Gain 

We consider the energy-to-peak gain response of linear time-invariant (LTI) systems. 
The specific application of interest is the computation of peak response bounds for robotic 
assembly of the International Space Station (ISS) subject to thruster firings. Mathematical 
modeling of the ISS robotic assembly results in high-order nonlinear models that contain linear 
structural dynamic effects and nonlinear robot arm dynamic effects. Since the operational 
scenario under consideration is static and the response of the robotic system is dominated by its 
local behavior, a linear representation is appropriate if the non-linear brakes-on behavior is 
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approximated by a linear equivalent. This can be accomplished by modeling the brakes-on 
dynamics by an equivalent linear stiffness. With this assumption, a linear time-invariant model 
from thruster inputs to joint torques can be obtained. Then, consider a linear time-invariant 
system, G(s), which has the state space form as follows, 

x = Ax + Bu 

y = Cx 

where u is the input pulse train. The objective is to estimate the peak pulse response of the above 
system without a time-domain simulation. This can be accomplished by making use of the 
energy-to-peak gain of the system. 

The peak transient response of a linear time-invariant system to bounded energy input 
signals can be computed from the energy-to-peak gain (Z^-to-L,*, gain) of the relevant transfer 
function [1]: 

r < P = su p IML (1) 

Hence, the energy-to-peak gain provides the worst-case peak response for the class of all unit- 
energy bounded signals. This gain can be computed from the following result [1,2] 

r,Hl cxc T 

where X is the controllability grammian of the system, which must satisfy the linear Lyapunov 
equation: 

AX + XA t +BB t = 0 

Hence, the energy-to-peak gain can be computed by solving a Lyapunov equation of order equal 
to the order n of the system. Efficient Bartels-Stewart algorithms can solve the Lyapunov 
equation requiring 0(n 2 ) floating point operations [3]. The symmetric solution can also be 
solved by use of the Kronecter product [4]. For very large-scale systems approximate solutions 
can be computed very efficiently using Krylov-subspace methods [5]. 

The above standard formulation of the energy-to-peak gain theory provides a tight upper 
bound as well as a method to construct the particular input signal for which the bound is reached. 
Hence, a straightforward application of this approach to a specific class of inputs, such as, 
thruster pulse trains can be very conservative. The conservatism is due to the fact that the 
standard formulation cannot, for example, distinguish between unit-norm bipolar and unit-norm 
unipolar pulses since both pulse trains have unit norm, resulting in the same estimate for the peak 
response. On the other hand, these signals generate very different time-domain signatures 
resulting in different peak responses. Thus, direct application of this method to practical 
problems involving specific classes of inputs would result in conservative estimates due the large 
class of input signals it allows. 

One approach to improve the peak response estimate is to restrict the allowable class of 
signals by preshaping the input signal so that it corresponds to a desired class of practical signals. 
From (1) it is apparent that the input signal u can be preshaped by a linear operator F resulting in 
a new input signal w, such that, u= Fw, as long as its induced 2-norm is less than one. In this 
case, the fictitious input w is any unit-energy bounded signal and F is selected so that Fw 
approximates the desired class of input signals. Hence, F can be used to shape or optimize the 
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resulting peak estimate. In this paper, a method to design such shaping filters is proposed to 
reduce the conservatism of the energy-to-peak gain for pulse train input signals. However, the 
methods can be generalized to any specified class of input signals. 

3. Shaping Filter Design 

The fundamental shaping filter design strategy is to exploit the frequency domain 
characterization of thruster pulse trains, specifically its power spectral density function, in order 
to design shaping filters that match the signal frequency content. The frequency domain 
representation of the input signals provides the ability to discriminate between similar norm 
signals. For the example presented in the previous section, the power spectral density (PSD) 
representation of unit norm bipolar and unipolar signals clearly allows for discrimination among 
these inputs. From the power spectrum of the input signal, filters can be designed such that the 
power spectrum of the filtered signal matches that of the input signal. 

Given a characteristic input or class of inputs in the time domain the corresponding filters 
can be designed by solving a frequency-based interpolation problem. First the Fast Fourier 
transform(s) (FFT's) of the input signals is calculated, and then a continuous-time filter, which 
corresponds to the FFT of the pulse train, is computed using any one of many standard 
interpolation algorithms. In our case we use the MATLAB function invfreqs to solve for the 
filter interpolation [6]. Assume the filter transfer function has the following form: 

, B(s) b n s n +b n l s n -'+...+b x s l +b Q 

r (5) = = 1 

A(s) a n s" + a n _ x s n + ... + a x s + a 0 

The coefficients of the numerator and denominator are determined by the minimal of the 
following cost function [6], 

n 

min^ wt(k) I F(k)A(k)~ B(k)\ 2 

a,b ^ 

*=1 

where wt(k) is a user specified weighting function. In our approach, only peak PSD response 
frequencies were selected and used in the solution of filter coefficients to avoid solving a large 
system equation. The magnitude at each peak response was selected to be the weighting 
function, wt(k). 

To guarantee existence of the Lyapunov solution [7], the filter has to be stable since the 
original system is stable. This can be implemented by removing the unstable poles of the filter. 
Once the filter has been formed, it is then normalized in order to achieve an induced 2-norm 
equal to 1 . The filter is then concatenated with the original system resulting in the augmented 
plant: 

G(s) = G(s)F(s ) 

where G(s) is the transfer function of the LTI plant and F(s) is the transfer function of the 
shaping filter. The augmented plant is then used to compute the energy-to-peak gain. 
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4. Test Case Results 


In this paper, the International Space Station assembly stage 4A is used to illustrate the 
the improvement in peak response estimate using the energy-to-peak gain formulation and 
shaping filters. The testcase represents the US Lab installation using the SRMS shown in Figure 
1. The testcase corresponds to the US Lab in the Pitch 180 configuration with the SRMS in 
brakes-on mode. The output variables of interest are the SRMS joint torques, shoulder yaw 
(SHY), shoulder pitch (SHP), elbow pitch (ELP), wrist pitch (WRP), wrist yaw (WRY), and 
wrist roll (WRR). The input signals are typical bipolar and unipolar thruster pulsetrains. The 2- 
pulse bipolar pulsetrain is shown in Figure 1 and consists of 20sec thruster firing time, 75sec 
delay between thruster firing, and the signal is sampled at 0.2sec. The PSD of the pulse train and 
the shaping filter are shown in Figure 3. The frequency response of the shaping filter is shown in 
Figure 5. The 5-pulse unipolar pulsetrain is shown in Figure 5 and consists of lsec thruster 
firing time, delays of 17.8, 1 1.4, 8.8, and 14sec between thruster firing, and the signal is sampled 
at 0.2sec. The PSD of the pulse train and the shaping filter are shown in Figure 6. The frequency 
response of the shaping filter is shown in Figure 7. 

The SRMS joint torque peak response is computed from simulation and compared to the 
estimates produced by the standard formulation (unweighted) and shaping filter formulation 
(weighted) of the energy-to-peak gain. The peak response is obtained for three pure single axis 
(roll, pitch, yaw) and a single multiaxis thruster pulsetrain testcases. The results for the bipolar 
inputs are shown in Table 1, while those for the unipolar input are given in Table 2. For the 
bipolar input testcase, Table 1, it is evident that the unweighted results substantially overestimate 
the peak response when compared to the simulation results. The overestimate is on average 
approximately 760% of the simulation peak response. The weighted results on the other hand 
reduce the overestimate by approximately 75%. For the unipolar testcase. Table 2, it is evident 
that the unweighted results also substantially overestimate the peak response when compared to 
the simulation results. The overestimate is on average approximately 860% of the simulation 
peak response. The weighted results reduce the overestimate by approximately 50%. Further, 
the overestimate in the weighted results can be attributed to errors in interpolating the frequency 
content of the input signal. 

5. Conclusion 

In this paper, shaping filters were proposed as a means to reduce the conservatism of 
analytical peak response estimates due to thruster pulse trains using the energy-to-peak gain 
formulation. The methodology was illustrated using two types of input signals for a realistic 
example involving Space Station robotic assembly. The results indicated that the use of shaping 
filters indeed does reduce the analytic peak response bounds substantially. Improved filter 
interpolation algorithms are expected to result in further improvements in the analytic estimates. 
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Figure 1. Stage 4A SRMS US Lab installation intermediate configurations 
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2-Pulse Bipolar Pulsetrain 
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Figure 2. 2-pulse Bipolar Signal 
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Table 1. 

Maximum Joint Torque 

for ISS 4A US Lab Pitch 180 due to Bi 

polar Input 

Inout 

OutDUt 

Sim 

Unwtd 

Wtd 

Unwtd/Sim 

Wtd/Sim 


SHY 

14 no 

76 07 

76 IQ 

S 44 

1 87 


SHP 

22.78 


41.64 

4.28 

1.83 

ROLL 

ELP 

16.69 

■ 

30.20 

4.20 

1.81 

WRP 

5.82 

I 

10.38 

4.30 

1.78 


WRY 

3.83 

26.58 

8.79 

6.95 



WRR 

3.15 

14.47 

; 5.43 

4.60 

1.73 


SHY 

4? Q4 

941 4? 

67 S7 

S 67 

1 46 


SHP 

31.41 

192.81 

55.98 

6.14 

1.78 

PITCH 

ELP 

22.85 


46.23 

10.96 

2.02 

WRP 




11.09 

2.07 


WRY 

16.43 

107.28 

23.84 

6.53 

1.45 


WRR 

12.08 

146.65 

23.47 

12.14 

1.94 



Averatw 


Table 2. Maximum Joint Torque for ISS 4A 


Outout I Sim I Unwtd 



due to Unipolar Input Signal 


Unwtd/Sim I Wtd/Sim 


Wtd 


3 6 iF-t.ni 
4.27E+01 
3.06E+01 
1.13E+01 
1.26E+01 
6.87E+00 


1 37F+07 
1.04E+02 
1 .32E+02 
L09E+02 
6.13E+01 
8.55E+01 


1 70F.+09. 
6.50E+01 
4.32E+01 
2.94E+01 
5.81E+01 
2.13E+01 


1 R7F.+07. 
L30E+02 
1.42E+02 
1.14E+02 
8.54E+01 
8.84E+01 
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